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DRIVING STANDING WAVES BY HEAT ADDITION 
By Perry L. Blackshear, Jr. 


SUMMARY 

Types of burner instability are enumerated and the role of standing 
waves in burners is discussed. The status of the problem of flame-driven 
standing waves is reviewed and a one-dimensional flow theory giving the 
mechanism whereby a flame drives or damps a standing wave is presented. 

In this theory, the reflection, transmission, and amplification of waves 
passing through a flame region were determined from the continuity and 
momentum equations. For the model considered, waves were found to pass 
through the flame front with their velocity amplitudes unaltered so long 
as the flame area remained unchanged. A change in flame area acted as a 
source of waves propagating simultaneously into the hot and cold gases 
on either side of the flame zone. 

1 - 2M 

0 

When the parameter — - — is unity (where Mq is average 

VTi/T 2 + 2Mq 

inlet Mach number and T 1 /T 2 the temperature ratio across the flame 
zone) , these waves are of equal magnitude and the driving criterion is 
that proposed by Rayleigh, namely, that for heat to drive a standing 
wave, the heat input should maximize with time at a place where the 
pressure in the standing wave varies and at a time when the pressure is 
near its local maximum value. 

1 - 2M o 

For — ■■ — >1/ the f lame-generated wave passing into the 

V T l/ T 2 + 2 Mo 

1 - 2Mq 

hot gas dominates and for — j = 1 < 1, the flame -generated wave 

W T 2 + 2M 0 

passing into the cold gas dominates. Phase requirements for driving and 
damping at these conditions are given. 

The role played by a flame in a standing-wave system was examined 
experimentally by measuring the ability of a flame to damp a standing 
wave. Factors investigated were: fuel-air ratio, inlet gas temperature, 

sound amplitude, inlet velocity, flame-holder position, and flame area 
as a function of time. 
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The percentage flame -area disturbance was found proportional to 
percentage velocity disturbance entering the flame zone. The phase lag 
of the area disturbance behind the velocity disturbance was found 
dependent upon flame-holder geometry and flame speed; in general, this 
lag increased as the flame speed decreased. 

The flame shape was strongly influenced by a radial component of 
the time -varying flow which appeared to follow a potential-flow velocity 
distribution in the neighborhood of the flame holder. This radial com- 
ponent of flow caused an initial growth of the disturbance in the flame 
front until a position was reached such that normal flame propagation 
caused these disturbances to undergo an apparent decay. 

The one -dimensional theory seems adequately to explain the experi- 
mental observations. 


INTRODUCTION 

In the field of jet-engine combustor research, burner instabilities 
are frequently encountered. A brief description of these instabilities 
is found in the names they have been given: in rockets, chugging and 

screaming; in afterburners, rumbling and squealing; in ram jets, rougji 
burning, resonance, pulsing, and whistling; even in Bunsen burners, 
occasional buzzing. 

In most cases a given combustor that exhibits one or more of these 
instabilities can be modified by a number of trial-and-error methods and 
the instabilities eliminated within the design range of operating condi- 
tions. For example, chugging rockets are quieted by increasing the pres- 
sure drop across the injectors; squealing in afterburners is subdued by 
increasing the velocity past the flame holders. 

The art of combustor design in its present state allows the 
designer, in most cases, to "cure" a specific burner of instabilities 
when they arise but can in no way guarantee that a new design will be 
free of them. Any of these instabilities can be extremely detrimental 
to engine performance. 

Fenn, Forney, and Garmon (reference l) have identified three types 
of instability found in ram- jet-type burners. They, along with a number 
of other investigators, have identified one type as resonance of the 
standing-wave or organ-pipe type. There is further evidence that this 
resonance can trigger other forms of instability that involve fuel- 
injection systems or supersonic diffusers. 

In 1859, Rijke (reference 2) discovered a method of driving standing 
waves by heat addition. Rayleigh (reference 3) and subsequently others 
have postulated an explanation of the Rijke tones (references 4 and 5), 
and have defined a criterion whereby flames (references 6 and 3) drive a 
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standing wave. In general, this criterion for driving a standing wave 
by the periodic addition of heat has been that heat be added at a place 
where the pressure in the standing wave varies and at a time when the 
pressure is a max imum . 

So far no mechanism whereby this heat addition can drive or damp an 
oscillation has been described in detail. An attempt to do this is made 
in the following theory of flame-driven standing waves derived at the 
NACA Lewis laboratory by employing one -dimensional flow equations to 
deduce the laws governing the transmission, reflection, and amplifica- 
tion of waves that approach a flame front. 

This point of view reveals the important parameters controlling the 
ability of a flame to drive or to damp an oscillation. The way these 
parameters vary in an actual flame is examined experimentally by meas- 
uring ability of the flame to damp a controlled imposed oscillation. 


Presuppose the existence of a standing wave and consider the region 
in a straight duct bounded on each end by a pressure loop at x = 0 and 
x = L, as shown in figure 1. (All symbols are defined in the appendix.) 
The flame holder is located at x^ at which point there is a flame 

region, assumed to be short compared with L. In the region O^x^Xf, 
the gas is at temperature T^_ and, in Xf <x < L, at temperature T 2 • 

Assume the wave equation applies to these regions considered separately. 
The boundary conditions connecting the two regions at x^ can now be 

obtained from the continuity and momentum equations. 

Consider the magnified view of the flame zone in figure 2 where the 
flame is shown anchored on a flame holder, and the zone 1 to 2 is consid- 
ered incompressible. The velocity entering plane 1 is taken to be 


THEORY OF FLAME -DRIVEN STANDING WAVES 


and leaving plane 2 to be 



( 1 ) 


The instantaneous flame area is 


r 


Af = Aav + a 


( 2 ) 
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so that 

u av A “ "f ( 3 ) 

where Uf is flame speed; U av , U ay Pp/P2> A^ are mean values; v-p 
v } and a are perturbations; and A is the cross-sectional area of 

u 

the duct. 


Since the region between 1 and 2 is incompressible, 

dv 

^1 = AfUf + ^ r 

and f 

Pi 3v f 

AU 2 = AfUf — + ^ J 

where is the volume of cold gas enveloped by the flame. 

From equations (l) to (4 ) , 


(4) 


U /Pi 
, av / 1 n 

V 0 a V. + a-T — 1 

2 1 ^av \ p 2 


(5) 


In order to obtain the pressure relation across the flame zone, 
consider the one- dimensional momentum equation for inviscid incompressible 
flow : 

p x . p 2 = P2 u 2 2 - Pl u x 2 + 1 dv 


= P2 U 2 2 ' Pl U l 2 + | §£ (Pl u l v f + p 2 U 2 V h) 


2 2 1 ( CP-f f 

W - Pl U l - 1 Ml sr - p 2 u 2 sr' + 


' iv 3plUl 
a \ Vf “5t 


+ V- 


*> 2 V 

h dt j 


( 6 ) 
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Where is the volume between planes 1 and 2 having density P2 and 

v h " ( x 2 " x l) A " v f 
But from equations (l) to (4) 

dv f c>V h au av 

Alt " Alt V 1 ' A av ' 

Inserting equations (l) to (5) and (7) in equation (6) and neg- 
lecting squares of perturbations gives the following expression for 
equation (6) : 


Pi - Pp - Pi ( — - 1) u 


av 


1 + 


2a 

*ai 


+ ll V f 


dp-,11 


1 U 1 


It" 


+ V 


Sp 2 U 2 

ll It" 


Now 


1 L 3P 1 U 1 SP 2 U 2 

A V Vf "St - + V h “5t 


can be shown from reference 7 to be approximately equal to 

dPi dp n U n \ ( Sp 

§7 


Ax n 


a „ dp i u i\ , 

u i —) + ^ 


2 a °2 D 2 
d's. ^ rjx 


These terms represent variously the gradients in pressure and momentum 
in the cold and hot portions of the flame zone which have been assumed 
short compared with wavelength and may be neglected. 

Outside the flame zone the wave equation applies. For the cold 
region (0<x<Xf) considered separately (reference 7), 


at 2 


= Cq 


dx^ 


the complete solution for which is 


h - Fk it 


C 1 


t + fk 


+ t 


This is for no flow, where \ is particle displacement, cj_ is speed 
of sound in 0<x<Xf> , and k is frequency. 
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The terms F and f are identified as right- and left -traveling 
waves, respectively, propagating at the speed of sound . For the case 

with flow, the right propagation velocity is c^ + U av and the left 
propagation velocity is Cj - XJ av ; | becomes 


!l=F 


k ( t Y1 + f k ( ^rr— + t) 

_ \ C 1 ^av U L v c l - ^av — 


Expanding 


C 1 + u av C 1 


= ~ H - Mq + M 0 


and 


“77 — = — ( 1 + M. + M„ 2 + . . 

C 1 - U av C 1 1 0 0 


where 


U. 

“Co 


av 


= inlet Mach number. 


Defining 


V = t + Mo | 

The general solutions may then be approximated by 



+ f 



Therefore, for moderate Mach numbers the succeeding arguments for 


apply equally well for t’ . 
used in x^<x<L.) 


(Since 




the same 


t 

t' 


is 


For small perturbations that satisfy the boundary condition at 
x = 0, namely, 0, the particular solutions may be written 


£]_ = B-l sin -t^ + Bj_ sin k^- + t^ => 2Bj_ cos kt sin ^-x (9) 
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where is amplitude. 


and 


The following equations are employed from reference 7 : 

2 S *1 
pi = - c i 2p i 77 

*1 

1 St 


( 10 ) 


(11) 


where and are perturbation pressures and velocities for this 

region and it is noted that for a right-traveling wave 

p 2 “ P 2°2 V 2 


and a left -traveling wave 


p i = - p i c i v i 


p 2 ’ - P 2°2 V 2 


r 


(12) 


J 


Similar relations are derived for x^<x<L such that the boundary con- 
dition at L is satisfied, again £ ^ = 0. 

The boundary conditions at x^ state that ' 


a / p l 

V 2 = Vl + V MpJ 




Let 


P 2 


a r 1 
U„„( — - 1 


( 5 ) 

( 8 ) 


Aav av \p2 


such that equation (5) becomes v. 3 v, + a* and equation (8) becomes 
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P_ = p, — - 1 U * + 2p^U a' 

2 !\P2 / av 1 av 

These equations are interpreted as follows. When a' is zero, 
waves pass through the flame zone with their velocity components unal- 
tered. The pressure components in this case are then determined by 
relations (12) . 

When a' is not zero, waves are generated at the flame front trav- 
eling into the hot and cold gases. The net change in velocity across the 
flame holder caused by these waves is equal to 


i ' - 


= V + V 

and a net change in pressure caused by these waves is equal to 

2pi u av a ' = - Pl 0 !^' " P2 c 2 a ? 
or, solving for a r * and a-^ ' , 



U. 


where 


av 


**0 "ST" 


= inlet Mach number and 


(13) 


(14) 


(15) 


a i' - - a ' 


P2 C 2 

Pl°I 


2Mq) 


1 + 


P 2 C 2 

p l c l' 





(16) 


where a r ' and a^ ' are the right- and left-traveling flame -generated 
waves. The ratio of magnitudes of these waves is 


1 - 2Mq 


1 - 2M, 


0 


a ' 


P 2 C 2 

Pl c l 


2Mq 



(17) 


2Mq 


* 


i 


If this parameter is greater than 1, the right-traveling wave will 
dominate } when it is less than 1, the left-traveling wave will dominate. 


2425 


2425 


2F 


NACA TN 2772 


9 


If the net effect is such that flame-generated waves are phased with the 
corresponding right- and left- traveling components of the standing wave, 
the flame will tend to drive the oscillation; when they are 180° out of 
phase, the flame will damp the oscillation. The net velocity v^ causes 

the disturbances in flame area. The phase lag of a' behind v^ where 
maximum driving or damping occurs is therefore of interest. 

For Vt as - 2Bk sin kt sin kx^/c-^, let 

a' = - H sin k(t-Cp) 

where H is a proportionality factor. 


But the right-traveling component of v.^ is 


Xf 


v^ = -Bk cos k\^— - tj = -Ek sin 




or for a ! to te in phase with 


cp = 2 rtx 
^ 2 k c 


k = 


2 * 


3 , x 
cp = 4 T + c 


Similarly, for a' to he in phase with v-^, 



Xf 

C 1 


Obviously, a' can supply right and left travelers simultaneously 
in phase with v^ and v^ at x^ = 0 or Xf = A/2 only. For any value 

of X£ and cp , the relative ability to drive for a given a 1 is 


D 


a 


' =c 


1 ~ R (cos kcp sin — 

1 + R V ci J 


kXf 

sin kcp cos — - 
C 1 
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where 


and for a ' cc v^_ 


1 

R 


1 - 2Mq 


P2 c 2 

P 1 C 1 


+ 2Mq 


1 - 2Mq 



D a 'covj 


kx~ 

sin ® a « aC 

C 1 


where negative values of D imply damping, and D is driving effective- 
ness. Maximum driving or damping occurs where 


cp = 



A plot of phase lags for maximum driving and maximum damping, respec- 
tively, cp + and cp_, in fractions of T as a function of Xf in 
fractions of is given in figure 3 for values of 


1 - 2Mq 

V T l/ T 2 + 2M 0 


and 2 


In addition, the envelope of constant D a ' , D a i oc v' , and a few rep- 
resentative values of D a t oo v' for a few arbitrary values of cp are 
plotted against x^ . 

In figure 4 relative driving effectiveness is plotted against cp 
for several values of x^. 

The driving criterion of Rayleigh, then, is for the case of 

1 ~ 2M o - x 
VVt 2 + 2Mo 

since for this condition, Cp + requires heat release to occur when pres- 
sure is a maximum. 

The experiments Rayleigh correlated with this criterion were for low 

M q and V T ]/ T 2 ~ 1 ‘ 


% 
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Other properties of the resonating system may be deduced. For cp 
where maximum driving or damping occurs, the period of oscillation will 
be 


or 


T = 


2Xf 2(L - Xf) 


Cl (l - Mq 2 ) 


1 - 


P-.U 

1 av 
P 2 C 2 


T « 2x f + 2(L - Xf) Vp 2 /pi 

c l(! - M q 2 ) Cl (l - M 0 2 Pl /p 2 ) 


For Cp other than cp + or cp_, the flame- generated waves will 

alter the wave form and the frequency. This frequency shift is dependent 
upon the departure of cp from cp + and upon the ratio of the amplitudes 

of the flame -generated waves to the existing standing waves. For any 
value of cp and for a ! DS.B, the frequency of resonance is 


1 "(irt sin_1 ir) sin k(cp + 

“ 2x1 2(L - x f ) 


cp) 


cid - ^ ) 


p,U 

1 av 

P 2 C 2 /J 


(where D is maximum possible effectiveness at a given x^) . Plots of 

relative driving effectiveness and relative frequency are shown in fig- 
ure 5 plotted against cp. These plots apply for any flame -holder posi- 
tion Xf and for any value of R. 


Interpre tat ion 

Before the experimental investigation is discussed it might he 
advisable to give a resume of the physical implications of the theory. 

For the model considered, waves can pass through the flame zone with 
their velocity components unaltered, providing that the flame area 
remains unchanged. For this case there is no driving or damping. If the 
flame area changes, the flame can drive or damp depending on the phase 
lag between the gas velocity and flame-area perturbations and the magni- 
tude of the area perturbation, as shown in figure 3. The driving 
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effectiveness given in figure 3 gives an index to the increase in wave 
amplitude each time it passes through the flame zone. If this increase in 
amplitude exactly compensates the loss per cycle to nonflame damping, the 
system will resonate. If this increase is greater or less than that 
required to compensate for losses the amplitude will increase or decrease, 
respectively, until flame "blow-out occurs or an equilibrium is reached. 

It was assumed that the flame zone was short compared with the 
quarter wavelength. In a duct it may he possible to excite the first 
several modes of vibration without violating this assumption. In this 
case the region of interest is the half -wavelength segment between the 
two pressure loops bracketing the flame zone, and figure 3 applies to this 
region. 

The mode of vibration the flame will select to drive will probably 
depend upon sound intensity and will probably be the one for which the 
ratio of flame amplification to duct damping is the greatest. 

Brief mention should be made concerning the point at which a flame- 
driven standing wave would give way to another type of oscillation. In 
brief, the standing wave should exist where flame amplification and sub- 
sequent attenuation per cycle is small. The point where this ceases to 
apply would be in the region where the attenuation at one end of the 
resonating system is such that the strength of the reflected wave equals 
the strength of the flame-generated wave alone. On the far side of this 
point, anticipated resonance is of a type having an acoustic length 
either two or four times the length from flame to reflecting end (as 
opposed to attenuating end) and a frequency dependent upon this acoustic 
length and upon the time lag between disturbance (velocity or pressure) 
and resulting flame -generated wave. 

EXPERIMENTAL INVESTIGATION 

It was felt that far more control could be exerted in examining 
experimentally the effect of a flame on a standing-wave system by meas- 
uring the ability of the flame to damp an Imposed standing wave than in 
examining spontaneous excitations. This was accomplished with the fol- 
lowing system: 

Apparatus . - A metered homogeneous mixture of propane and air was 
fed to a 1- inch- diameter burner as shown in figure 6, through a variable- 
area pintle-type nozzle. Standing waves in this burner were driven by a 
piston speaker through a high- impedance source tube (reference 8) at the 
base of the burner. Sound level was measured in the burner base by two * 

flush-mounted button microphones. The flame holders employed (fig* 7) 
could be positioned in a test section where desired. Auxiliary equipment 
included a hot-wire anemometer, a stroboscope, a camera, and a photocell. 
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Properties of system . - The Reynolds number in the test section was 
about 500. The velocity profile with sound appeared to be a superposition 
of a potential-flow profile for the time-varying flow upon the steady 
laminar flow (fig. 8) with no flame holder present. The time- varying 
velocity at points across the test section is indicated by the numbered 
inserts, showing that near the walls flow reversal occurred. Placing a 
flame holder in the test section gave velocity profiles shown in fig- 
ure 9. A comparison of the time-varying component with a computed (ref- 
erence 9) two-dimensional potential flow in figure 10 shows that fair 
agreement exists except for a boundary layer. This boundary-layer thick- 
ness seemed to decrease with an increase in Reynolds number of the steady 
flow. 


A variety of mufflers were examined and one was selected having a 
fairly constant attenuation for the first four modes of vibration. 

As the pressure drop across the inlet orifice increased, the random 
noise increased while the damping coefficient of the burner decreased 
until a pressure drop of 1.5 inches of mercury was reached. Above this 
value, the noise continued to increase but the damping coefficient 
remained substantially constant. 

A sample of the spontaneously excited frequencies found in the ini- 
tial exploratory work is shown in figure 11. Here a 64-mesh-per-inch 
screen was used as a flame holder. The exhaust end of the tube was 
alternately open or masked by a plate containing a l/l6- inch-diameter 
orifice. With this arrangement, that mode of vibration was excited such 
that the flame holder was upstream of a pressure loop, with the tube 
acting as either a closed-closed (X/2) or closed-open (X/4) resonator. 

Procedure. - Two indexes of damping coefficient were employed (ref- 
erence 10) : the width of resonant peak, or the voltage required to drive 

the speaker for a constant sound intensity in the burner. The role of 
the flame in damping the imposed oscillation was investigated by consider- 
ing the following variables: fuel-air ratio, inlet temperature, flame- 

holder position, sound amplitude, and inlet velocity. 

Effect of fuel-air ratio on damping. - The effect of fuel-air ratio 
on damping is illustrated in figures 12(a) to 12(c), for the first four 
modes of vibration with all three flame holders. The curves seem to be, 
within limits, symmetrical about the fuel-air ratio of max im u m flame 
speed, 0.07 (reference 11) . That this would indicate a strong dependence 
on flame speed is not immediately evident since flame-speed studies (ref- 
erences 11 and 12) show flame speeds as asymmetrical about a fuel-air 
ratio of 0.07. 

However, these references do not agree on flame speeds for rich mix- 
tures; reference 12 indicates higher and reference 11 lower values than 
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those required for symmetry. A few check readings of photographs of 
flames 'burning on the one -hole flame holder within the tube showed symme- 
try about a fuel-air ratio of 0.07 and agreed moderately well with the 
reference lean values. 


Tn figure 12(d) the voltage to speaker and the observed frequency at 
resonance for the third mode of vibration are plotted against fuel-air 
ratio. For the three- and five-hole flame holders, these data show good 
agreement with figure 5 if fuel-air ratio is taken as an index of cp 
such that cp increases as flame speed decreases . 


This is interpreted to mean that the phase lag of the flame-area 
change behind the velocity perturbation has a strong flame-speed depend- 
ence. A comparison of figures 12(a), (b) , and (c) indicates a strong 
dependence on flame-holder geometry. 

Data shown in figure 12 were taken in a steel test section. Some of 
the heat transferred from flame to wall downstream of the flame holder 
was conducted upstream and thus preheated the combustible mixture. There- 
fore, when these figures are compared with subsequent curves taken in 
glass tubes, this preheating must be considered. 


Effect of inlet temperature on damping . - The effect of fuel-air 
ratio on damping is shown in figure 13 for various gas temperatures. The 
fuel-air ratio at which damping peaked for each temperature was such that 
flame speed remained the same. If a/^av is assumed proportional to 


v^/Q , and sound intensity is assumed constant 



then 


a’ 



CD v-l 



A slight increase in the driving or the damping as temperature increases 
would therefore be expected for a constant mass flow, as indicated in 
figure 12. 

Effect of flame-holder position on damping . - The effect of flame- 
holder position y on damping coefficients for the third mode of vibra- 
tion is shown in figure 14. It has been shown from theoretical considera- 
tions that the damping would be expected to depend upon a/A av , hence 

v 1 , and upon the phase lag cp. For the third mode a marked peak is 

obtained in the curve of damping against fuel-air ratio curve for each 
position except at y = 0. At y = 0, the right-traveling velocity 
wave will not be reflected, hence will not be able to drive or damp, and 
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the left-traveling flame-generated wave will he weaker because of a 
diminished momentum pressure drop. In the range from y = 2 to 

y = inches , the magnitude of flame damping increased and the fuel-air 

ratio where maximum damping occurred shifted toward the fuel-air ratio 
of maximum flame speed. This effect would imply from theory that phase 
lag diminishes for the third mode as flame speed increases , as is demon- 
strated subsequently. Maximum damping would be expected to occur at 
y = }\/Q * 10 inches for this mode. The spontaneous excitation of a 

higher harmonic at y >8^ inches precluded verification of this maximum 


Effect of sound amplitude on damping . - Over the range of amplitude 
examined (120 to 160 db; , there was no perceptible change in the damping 
characteristics of the flame, as long as the flame remained seated. This 
would indicate that the flame component of damping varied directly as the 
disturbance velocity, that is, a/A av co /U av ' Furthermore , it would 

imply that amplitude did not alter the phase lag. Too large an amplitude 
would cause the flame to lift or flash back. Tests run at amplitudes 
that were not large compared with the background noise were impossible to 
interpret with the equipment employed. A number of strange effects were 
noted. For example, with a combustible mixture flowing through the tube, 
a flame ignited at the open end would bum upstream to within an inch of 
the flame holder and hover there, exciting spontaneously the first mode 
of vibration. Imposing the second- or third-mode frequency at a high 
amplitude would cause the flame to drop down to the flame holder, where 
it would seat and burn quietly after removal of the imposed vibration. 

Effect of inlet velocity on damping . - Damping against fuel-air 
ratio for different velocities is shown in figure 15 for the three -hole 
flame holder. The magnitude of damping is seen to vary directly as 
U av , whereas the curve shape remains essentially unaltered. This is 

interpreted as meaning the phase lag is essentially unaltered by flow 
velocity. 


would remain constant as U av is varied. The apparent increase in a' 
with an increase in U av is believed due in part to the Reynolds number 

effect of the steady flow upon the time-varying component of flow men- 
tioned earlier. The term a/A^, then, is more nearly proportional to 


point . 


If a/A av were proportional to v^ /u av , a 




where n in this instance is slightly less than unity. 
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Flame area as a function of time . - Figure 16 shows sample photo- 
graphs of periodically varying flames at different phases of velocity 
for the first and third modes of vibration. Here the radial distortion 
of the flame at the flame holder due to the radial component of the time- 
varying potential flow is apparent in the photographs of the third mode 
of vibration. A particle of gas near the flame-holder rim is given a 
radial momentum when it passes through the flame holder during outflow 
which it never loses in the ensuing inflow part of the cycle. This 
causes an apparent growth in the disturbance when the flame-propagation 
rate is less than this outflow rate. Because of normal flame propaga- 
tion, the flame eventually becomes sharply convex toward the burned gas 
at these outermost loci and then propagates at a rate that is locally 
greater than normal laminar flame speed. This gives the appearance of a 
flame-front disturbance, the amplitude of which is reduced from this 
height to the tip. The rapid transition to cycloidal cusps in the flame 
front appeared to be due solely to normal flame propagation in lean and 
rich flames, and not to flame-front instability (references 14 and 15) . 
That the disturbances did not grow as they have been shown to do in tur- 
bulent flames (references 13 and 16) would indicate that a/A. av , as 

observed here, is a lower limit. 


The first-mode photographs gave evidence of a different flow field. 
When the amplitude became sufficiently great, the flow appeared composed 
of successive ring vortices detached from the flame holder and disturbed 
the flame in so doing. The transition from this vortex flow to the 
potential flow of the third mode seemed to occur in the neighborhood of 
200 cycles per second for the conditions studied. 

Photographic sequences for the three-hole flame holder were taken 
over a range of fuel-air ratios bracketing the fuel-air ratio where 
damping maximized. 


Damping for these sequences and the frequency variations are shown 
in figure 17. Figure 18 shows relative instantaneous flame areas of 
these sequences plotted against degrees behind maximum velocity, giving 
a measurable lag of flame area, hence heat release, behind velocity. 
Maximum damping occurred at a fuel-air ratio of 0.0775. The phase lag 
for this fuel-air ratio is in the neighborhood of 130°. This is the 
phase lag which would be predicted from the theory for 


1 - 2M o 

V T l/ T 2 + 2 Mo 


2.2 


and for = 37^/16, which seems to fit the conditions imposed for this 

sequence. In addition, the phase lag increased as flame speed decreased 
so that the frequency shift is as predicted in figure 5. 
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For the three conditions studied photographically, the phase lag 
decreased with increasing lip for the third mode, and increased, then 

decreased with increasing for the first mode. No explanation for 

this devious coupling between flame area and velocity is immediately 
obvious . 


Summary of experimental results . - From the preceding examinations, 
a fairly clear picture of the relation between an inlet velocity pertur- 
bation and the resultant flame-area perturbation may be obtained. For 
the third mode and three-hole flame holder, for example, the magnitude 
of the area disturbance a/A^ is proportional to VjJV av n where n 

is slightly less than unity because of a boundary- layer effect. The 
phase lag of a ]_/A ay behind v-|_ depends only upon flame speed U f and 

not upon v^; this phase lag increases with diminishing flame speed. 

With these properties in mind, it is seen that the agreement with theory 
is quite good. Disparity exists between 'theory derived herein and the 

1 - 2Mq 

Rayleigh criterion when _ — ^ 1. The shortcomings of the 


Rayleigh criterion can be verified by experiment: (l) by obtaining a 

measurable phase lag where maximum damping occurs 


W T 2 + 2M 0 

: ( ; . 

s; (2) by showing a 
change in the apparent phase lag with x^; or (3) by a measured existence 

of driving or damping at a velocity loop. The experimental verification 
of these three phenomena is discussed in the following paragraph. 

In figure 18 the measured phase lag corresponding to maximum drive 
is seen to be approximately 130°, which would correspond to 

16' 


VT;l/T 2 


*2.2 and x„ « 


2Mq 


According to the Rayleigh criterion, 


this phase lag could be only 90 u in the region 0>x 4 ,<^i. These data 

— i — 4 

are not wholly reliable since the time-varying velocity is deduced from 
the flame slope at the base and would not account for a lag in the jet 
boundary. Subsequent checks made with a photocell and a hot-wire ane- 
mometer coupled to a phase meter verified the lags obtained photo- 
graphically. 

A plot of volts to speaker as an index of damping against fuel-air 
ratio as an index of phase lag for different flame-holder positions 
(fig. 14) shows a shift in the phase lag where damping maximizes that is 

1 - 2M 0 

consistent with the theory for the case of — — — > l and could 


Vt i/t 2 


2Mq 
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not be explained by the Rayleigh criterion. Finally, cases of the spon- 
taneous excitation of the fundamental mode of the duct were recorded with 
the flame seated at the open end of the duct. 


CONCLUSIONS 

The one -dimensional analysis seems an adequate explanation of the 
mechanism whereby a flame drives or damps a standing wave. This mech- 
anism is as follows : Waves pass through the flame front with their 

velocity amplitude unaltered. When the net velocity causes a change in 
flame area, waves are propagated simultaneously into the hot and cold 
gases. When these flame -generated waves are phased with the existing 
components of the standing waves, they drive the standing wave; when 
opposed, they damp. 


When the parameter 


1 " 2M 0 

— — is unity (where M n 

V T l/ T 2 + 2M 0 


is average 


inlet Mach number and T^/T^ the temperature ratio across the flame 

zone), these waves are of equal magnitude and the driving criterion is 
that proposed by Rayleigh, namely, that for heat to drive a standing 
wave, the heat input should maximize with time at a place where the pres- 
sure in the standing wave varies and at a time when the pressure is near 
its local maximum value . 


The strength of the wave depends largely on the relative magnitude 
of the disturbance of the flame area a/A^ and the phase relation 

between a/A ay and perturbation velocity v 1 . The magnitude of a/A av 
was found to depend upon v 1 and average flow velocity in the cold gas 
U ay and in most cases varied directly as v i/u av » The existence of 

radial flow in the time -varying component of the flow through the flame 
holders caused a greater perturbation than would be anticipated in one- 
dimensional flow. The phasing of a/A with v 1 was found to depend 

primarily on flame speed and, for the conditions studied, appeared 
independent of the magnitude of v^ /U av • In most cases, the phase lag 

of a/A behind v, /U increased as the flame speed decreased. 

' av 1' av 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 13, 1952 


2425 


2425 


NACA TN 2772 


19 


A 


Af 

■^av 

a 

a' 

V 

a r * 

B 

c 

D 

E 

F,f 

g 

k 

L 


“o 

n 

P 

p a. 

P 


APPENDIX - SYMBOLS 


The following symbols are used in this report : 
duct cross sectional area, sq ft 
instantaneous flame area, sq. ft 
average flame area, sq ft 
" Aav 

ft/sec 


left-traveling component of flame- generated wave, ft/sec 

right-traveling component of flame -generated wave, ft/sec 

amplitude of one component of a standing wave, ft 

velocity of sound, ft/sec 

driving effectiveness, dimensionless 

voltage required to drive speaker, volts 

arbitrary functions 

gravitational constant, ft/sec^ 

frequency = 2 «/t, radians/sec 

distance between pressure loops = X/2, ft 


= inlet Mach number 

C 1 


integer 1, 2, . . . 

instantaneous static pressure = P ay + p, lb/sq ft 

average static pressure, lb/sq ft 

static pressure of standing wave, lb/sq ft 
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R 


T 

t 

U = 


Uf 

V 

v c 

V f 


x 


y 

r 

A 



p 

T 

cp 


V T l/ T 2 + 2 Mo 

1 - 2Mq 

temperature 
time, sec 
U + v, ft/sec 

average flow velocity in cold gas, ft/sec 
normal laminar flame speed, ft/sec 
volume 

volume within flame zone at temperature 1, cu ft 
volume of cold gas enveloped hy flame, cu ft 
volume within flame zone at temperature 2, cu ft 
perturbation velocity, ft/sec 

distance in direction of flow from pressure loop, ft or in. 
flame -holder position measured from pressure loop, ft or in. 
flame-holder position measured from exhaust end, ft or in. 
ratio of specific heats, c p/ c v 
wave length, ft 
TC-p ft 

Tc 2 , ft 

displacement of particle in wave, ft 
average density, slugs/cu ft 
period, sec 

phase lag of a behind v^, sec 
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phase lag for maximum drive, sec 
cp phase lag for maximum damping, sec 

Subscripts : 

1 inlet to flame zone 

2 outlet to flame zone 
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Figure 1. - Segment of duct tietween two pressure loops. 



Driving effectiveness, D Phase lag 
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(a) Phase lag where maximum driving or 


damping 


1 - 2M 0 

occurs; - — ■— ^ — 

VV T 2 + 2M 0 


1 

2’ 



(b) Envelopes of driving effectiveness for 
a’cov^ and a' = constant; 



Flame-holder position measured from pressure 
loop, Xf 

(c) Driving effectiveness for representa- 



tive values of cp; _ H =. JL. 

-JT1/T2 + 2M 0 2 

(Numbers beside curves are values of <P.) 

Figure 3. - Phase-lag requirements for maxi- 
mum effectiveness. 
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a'oovj and a' = constant; 



Flame-holder position measured from pressure 

lOOp, Xf 

(f) Driving effectiveness for representative 


1 - 2M n 

values of cp; x * 1. 

V T l/ T 2 + 2M 0 

beside curves are values of ) 


(Numbers 


Figure 3. - Continued. Phase-lag requirements 
for maximum effectiveness. 
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(g) Phase lag where maximum driving or 


1 - 2M n 

damping occurs; X — = 2. 

V / T 1 /T 2 + 2M 0 



Flame-holder position measured from pressure 
loop, Xf 


(l) Driving effectiveness for representative 


values 

1 - 2M 

of q>; — = 2. 

/V T 2 + 2Mq 

(Numbers 

- 

beside 

curves are values of cp. ) 




Figure 3. - Concluded. Phase-lag requirements 
for maximum effectiveness. 
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, x 1 - 2M 0 

(c) 

yr i/T 2 + 2 Mq 

Figure 4. - Relative effectiveness at given 
x f for variations in op. (Numbers beside 
curves are values of x^.) 
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(b) Relative frequency at any x^ for variations in cp. 


Figure 5. - Relative driving effectiveness and frequency for any x^, 
against cp relative to for the chosen. 
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Figure 6. - Test apparatus. 
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Figure 7 . - Flame holders . 
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1.4 


1.3 


1.2 


1.1 


1.0 




© Velocity trace with 

steady and time- 


varying flow 



.2 .1 0 .1 
Distance from center line of duct, in. 

Figure 8. - Velocity profiles in test section with no flame holder present. Average flow 
velocity U av , 0.72 feet per second; sound level, 132 decibels at 22 cycles per second. 
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Figure 9. - Velocity profiles measured l/l6 inch downstream of one- 
hole flame holder placed 4 inches from exhaust at 55 cycles per 
second. 


2425 


Velocity maxima, ft/sec 


5F 


MCA TN 2772 


33 



Figure 10. - Comparison of measured time-varying flow through flame holder 
(see fig. 9) with two-dimensional potential flow. 
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Figure 11. - Effect of flame-holder position on spontaneous sound for exhaust end of burner acoustically 
open (V 4 ) or acoustically closed (\/2) • Over-all burner length, 58 inches; average flow velocity U av , 
0.875 feet per second; fuel-air ratio, 0.06; flame holder, 64-mesh screen. 
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Figure 12. - Voltage to speaker required to drive 
■burner at 123.2-decibel sound level at resonance 
for first four modes of vibration. Flame -holder 
position measured from ekhaust end y, 4 inches; 
average flow velocity U av , 1.51 feet per second. 



(b) Three-hole flame holder. 


Figure 12. - Continued. Voltage to speaker re- 
quired to drive burner at 123.2-decibel sound 
level at resonance for first four modes of vi- 
bration. Flame -holder position measured from 
exhaust end y, 4 inches; average flow velocity 
U ay , 1.51 feet per second. 
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Fuel-air ratio, f/a 


(c) Five-hole flame holder. 

Figure 12. - Continued. Voltage to speaker required to drive 
burner at 123.2-decibel sound level at resonance for first 
four modes of vibration. Flame-holder position measured 
from exhaust end y, 4 inches; average flame velocity U av , 
1.51 feet per second. 
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Figure 12. - Concluded. Voltage to speaker required to 
drive burner at 123.2-decibel sound level at resonance 
for first four modes of vibration. Flame-holder posi- 
tion measured from exhaust end y, 4 inches; average 
flow velocity U av , 1.51 feet per second. 


Voltage to speaker , E Voltage to speaker 
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(a) Third mode. 
44 t- 



.05 .06 .07 .08 .09 

Fuel-air ratio, f/a 

(b) Fourth mode. 


Figure 13. - Effect of inlet temperature on damp- 
ing for three-hole flame holder in glass tube. 
Flame-holder position measured from exhaust end 
y, 4 inches; average flow velocity U av , 

1.51 feet per second; sound level, 123 decibels. 
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Figure 14. - Effect of flame -holder position on 
damping for three-hole flame holder. Third 
mode . 



Figure 15. - Effect of inlet velocity 
on damping for three -hole flame 
holder in glass tube. Third mode; 
flame-holder position measured from 
exhaust end y, 4 inches; sound 
level, 123 decibels. 
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(a) Frequency, 329 cycles per second; voltage, 23.8 volts; velocity, 1.378 feet per second; 

fuel-air ratio; 0.080. 

Figure 16. - Sequence of flame shapes for one cycle at intervals of 45°. 

clockwise. 


Time increases 
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(b) Frb,««or, 3U cyolsa per 25. ,3 volt.; v.loolty, 1.37S feet per .pcoM; 

v ' fuel-air ratio, 0.090. 


Figure 16, 


- Continued. Sequence of flame shapes for one cycle at intervals of 45°. Time 

increases clockwise. 
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(c) Frequency, 318 cycles per second; voltage, 32.0 volts; velocity, 1.378 feet per second; 

fuel-air ratio, 0.0775; run 2 (fig. 17). 

Sequence of flame shapes for one cycle at intervals of 45°. Time 
increases clockwise. 


Figure 16. - Continued. 
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(d) Frequency, 58.5 cycles per second; voltage, 17.3 volts; velocity, 1.378 feet per second 

fuel-air ratio, 0.080. 

Figure 16. - Continued. Sequence of flame shapes for one cycle at intervals of 45°. 

increases clockwise. 


Time 
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(e) Frequency, 55 cycles per second; voltage, 12.9 volts; velocity, 1.378 feet per second; 

fuel-air ratio, 0.090. 

Sequence of flame shapes for one cycle at intervals of 45°. Time 
increases clockwise. 


Figure 16. - Concluded. 
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Fuel-air ratio, f/a 


Figure 17 . - Damping and frequency distribution for 
three -hole flame holder. Third mode. (Numbered 
inserts correspond to points where phase lag was 
determined from photographic sequences.) 
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Figure 18. - Flame area as function of time for three-hole flame holder. Third mode. 
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